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1. Motivation and Objective 

The Yakhot-Orszag renormalization group is used to analyze the pressure gradient- 
velocity correlation and return to isotropy terms in the Reynolds stress transport 
equations. The perturbation series for the relevant correlations, evaluated to lowest 
order in the e-expansion of the Yakhot-Orszag theory, are infinite series in tensor 
product powers of the mean velocity gradient and its transpose. Formal lowest 
order Pad£ approximations to the sums of these series produce a rapid pressure 
strain model of the form proposed by Launder, Reece, and Rodi, and a return to 
isotropy model of the form proposed by Rotta. In both cases, the model constants 
are computed theoretically. The predicted Reynolds stress ratios in simple shear 
flows are evaluated and compared with experimental data. The possibility is dis- 
cussed of deriving higher order nonlinear models by approximating the sums more 
accurately. 

The Yakhot-Orszag renormalization group provides a systematic procedure for de- 
riving turbulence models. Typical applications have included theoretical derivation 
of the universal constants of isotropic turbulence theory, such as the Kolmogorov 
constant, and derivation of two equation models, again with theoretically com- 
puted constants and low Reynolds number forms of the equations. Recent work 
has applied this formalism to Reynolds stress modeling, previously in the form of a 
nonlinear eddy viscosity representation of the Reynolds stresses, which can be used 
to model the simplest normal stress effects. The present work attempts to apply 
the Yakhot-Orszag formalism to Reynolds stress transport modeling. 

2. Work Accomplished 

The modelling of the pressure gradient- velocity correlation and return to isotropy 
term in the Reynolds stress transport equation remains an area of active research. 1 ,2,3 

Models will be derived here using the Yakhot-Orszag renormalization group 2 * 4 par- 

tially along the lines of our previous work 5 . The result is a model for the rapid 

pressure-strain term of the form proposed by Launder, Reece and Rodi 6 (LRR) and 

a model for return to isotropy of the form proposed by Rotta 7 with theoretically 

computed constants in good agreement with accepted values. As is usual in inves- 
tigations of this sort, the priority of Yoshizawa in deriving a pressure strain model 
analytically 8 must be noted. 

The analysis requires some new ideas in renormalization group theory recently 
introduced by Yakhot et al 9 . As Yakhot et al 9 emphasize, the application of the 
renormalization group mode elimination formalism to shear flow creates a double 
perturbation series in powers of e, the parameter of the isotropic theory, and in 
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powers of a dimensionless strain rate, 77 = SK/e , where K denotes the turbulence 
kinetic energy, e denotes the dissipation rate, and 5 is a measure of the mean 
strain: in Ref. 9, S 2 = (^ BL + The present analysis also leads to double 

expansions of this type, with the powers S n replaced by tensors homogeneous of 
degree n in the mean velocity gradient Vf7 and its transpose. It will be convenient to 
retain the terminology of Ref. 9 and call this expansion the 77-expansion; when the 
distinction is pertinent, the expansion of Ref. 9 will be called a scalar 77-expansion. 

The heuristic program of evaluating all scalar amplitudes to lowest order in e 
has proven successful in the past: apparently, the e-expansion is an asymptotic 
series with a sum given very nearly by its first term 10 . Unfortunately, there is no 
analogous basis for truncating the 77-expansion. There are fundamental reasons for 
this distinction between these expansions. The present 77-expansion is tensorial: 
successively higher order terms do not introduce merely numerical corrections, but 
increasingly complex asymmetries into the theory. Truncation therefore imposes 
a possibly inappropriate symmetry or other constraint on the model. Thus, in 
Ref. 5 the 77-expansion for the Reynolds stress r was truncated at second order as 
suggested by previous work of Yoshizawa 11 and Speziale 12 . Although this type of 
modelling permits unequal normal stresses in a simple shear flow, it is not maximally 
asymmetric: for example, in a flow with mean velocity components Ui(x i,x 2 )» a 
cubic model including a term r ~ VU 2 VU T + VUVU T2 would permit nonzero r 2 3 
in the presence of vanishing dU-i / dx$ and dUz / dx% , an effect which cannot be ruled 
out in advance. 

Although generalizations 13 of the Cayley-Hamilton Theorem limit the number 
of independent tensors S (n) , anisotropy and asymmetry cannot exist at all without 
some terms of higher order in 77; indeed, truncation at lowest order in 77 just produces 
a theory of isotropic turbulence. But the series truncated at any higher order can 
be unsatisfactory in flow regions in which some components {^)K/e are large. In 
such regions, the truncated series is dominated by its highest order terms. For the 
quadratic stress models of Refs. 5, 11, 12, this domination can produce negative 
normal stresses in the buffer layers of wall bounded flows. Increasing the order of 
truncation obviously exacerbates this problem. 

It follows that finite truncation of the 77-expansion is theoretically unsatisfactory. 
Yakhot et al 9 therefore propose that this expansion must be summed, even if only 
approximately, and have suggested a prototype summation in a different context. 
It should be noted that the same issues arise naturally in Yoshizawa’s formalism, 
which also generates infinite series in the mean velocity gradients (and in other 
quantities as well) for correlations of interest in turbulence modeling. Yoshizawa 
has concluded independently that summation of this series is essential and has also 
derived a Reynolds stress transport model by introducing such summations 14 . 

In this paper, the perturbation series which the Yakhot-Orszag renormalization 
group generates for the correlation 


IJij — 
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is summed by a low order Pade approximation. Coefficients 'are evaluated to lowest 
order in the e expansion, but the summation includes effects of all orders in tj. le 
result is essentially identical to the “model 1” proposed by Launder, Reece and 
Rodi 6 . An entirely analogous treatment of return to isotropy yields a model of the 
form proposed by Rotta 7 . Combining these models leads to a preliminary Reyno s 
stress transport model. The problem of closing the Reynolds stress diffusion terms 
is addressed. This problem also leads to an infinite sum. . 

While it is encouraging that renormalization group methods can be used to deny 
familiar models, the goal of this investigation is not limited to providing theoreti- 
cal justification for the LRR and Rotta models, which although widely applied are 
nevertheless deficient in several well-documented respects • ’ • Instead, renorm 
ization group methods together with approximate summation of the ^-expansion 
can be used to derive higher order and nonlinear corrections to these models in a 
systematic fashion. Explicit development of such models is left to future investiga- 
tions. 

2.1 Analysis of the Pressure Gradient- Velocity Correlation 

The analysis will foUow Yakhot and Orszag’s derivation of turbulence transport 
models by renormalization group methods 4 . The equation for velocity products is 


d 

m u ' Uj + 


( duj duA ( dp dp_\ 

U > { Ui dT p + Uj d^J " " r&r; +Uj 8xJ 


+i/o V 2 UiUj — 2vojj-^ 


dui duj 


dip dx p 


( 2 ) 


where i/ 0 denotes the kinematic viscosity. The product -( mdp/dxj + Ujdp/dxi ) 
on the right side of Eq. (2) will become the correlation defined by Eq. (1) 
following elimination of all fluctuating modes. 

Thus, the perturbation series will be written as 

II = To + Ti -I 

where T n is of order n in u< and all amplitudes are evaluated to lowest order in e. 
To lowest order in e and SK/e 




( 3 ) 


in agreement with the analysis of Crow. 16 

At the next order in SK/e , in the high Reynolds number limit, 
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where (0) denotes deviatoric part and ( ij ) denotes index interchange in the preced- 
ing term. 

The next order will produce a term T 3 containing cubic products of velocities u K . 
In view of the form of the LRR model, it is reasonable to ask whether a term with 
only one gradient, proportional in the high Reynolds number limit to tVU might 
occur at this order. Such terms do occur, but they cancel. Evaluation of T 3 proves 
to require expansions of the projection operators to second order, leading instead 
to terms S (3) homogeneous of degree three in the mean velocity gradient and its 
transpose. In general, the term T n of order n has the form S^(K/e) n . As noted 
in the Introduction, it will be imperative to include effects of all orders in SK/e 
in the model, but because the terms T n involve ever higher order derivatives of 
the transverse projection operators, they do not have an obvious law of formation. 
Therefore, an exact summation does not appear feasible. 

A simple approximate summation is obtained by introducing the perturbation 
series 5 for uiUj^ in the form 


v 




-Wuj {0) + Y, S (n) (*7<0 n 

n> 1 


and dropping the quadratic terms. The resulting model, 
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with 


-7619 


Cfi = ^ = 0952 


( 6 ) 


agrees with the perturbation series (3) and (4) to terms of order S (3) . However, 
unlike the explicit quadratic model which results from simply dropping the 0(S (3) ) 
terms, this model includes effects of all order in SK/e. The consequences of this 
fact will be discussed later. This type of summation has also been applied by 
Yoshizawa 14 . Eqs. (5) and (6) can be compared with Launder, Reece and Rodi’s 
“model 1”, Eq. (5) with the empirically adjusted constants 


Cf i = .7636 Cfi = .1091 (7) 

In this model, the constants Cf X and Cf 2 were not chosen independently; instead, 
to insure some consistency conditions introduced by Rotta 7 , Launder, Reece, and 
Rodi set 6 


C 


fi 


C*2 + 8 
11 


Cf 2 = 


8C 2 -2 
11 
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where only the constant C2 is arbitrary. By eliminating C2 between these equations, 
there results 


8CV1 -Cr 2=6 ( 8 ) 

which is also satisfied by the choice of constants in Eq. (6). The LRR model 
corresponds to the choice C 2 = .4; Eqs. (5) and (6) correspond instead to the 

choice C2 =8/21 ~ .36. 

The approximate summation used to derive Eq. (5) can be systematically gener- 
alized to generate an infinite number of models for Ilij. For example, suppose that 
the perturbation series for r is introduced into the cubic terms in the perturbation 
series instead of in the quadratic terms as above. This substitution will produce a 
model which can be written symbolically in the form 

iT-SW+S^+^'+SW) 

where rS^' denotes a sum of matrix products in all possible orders of r and terms 
S(‘). The requirement that the original series agree to order S (4) with the approxi- 
mation when r is replaced by its perturbation series determines this approximation 

uniquely. 

2.2 The Return to Isotropy Model 

The analytical description of return to isotropy is no less controversial than the 
modeling of the fast pressure strain term 3 . In the usual approach to turbulence 
modeling, in which correlations generated by Reynolds averaging are closed phe- 
nomenologically, this process is considered to result partly from the pressure cor- 
relation through a “slow” term independent of the mean flow, and partly from the 
deviatoric part of the dissipative correlation From this viewpoint, the 

analysis in Sect. I is incomplete because it discloses only a term proportional to the 
mean velocity gradient, but no slow term. The return to isotropy will be derived 
here by renormalization group methods following a suggestion of Yakhot . 

From the renormalization group viewpoint, it is natural to investigate the return 
to isotropy, even independently of the stress transport equation, by writing the 
perturbation series for 

Ui ln + I “ i( * " + W) < 9 > 

This perturbation series differs from the perturbation series for the Reynolds stresses 
previously reported 5 only in the occurence of an additional factor -iu in all fre- 
quency integrals. 

The analysis is straightforward. Only the deviatoric terms require attention be- 
cause the part of the correlation proporational to 6ij contributes to the transport 
equation for K which has been analyzed by Yakhot and Smith 15 . The lowest order 
deviator appears at first order in 77; to lowest order in e 
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where 


Ti 



dUj 

dxj 


+ 


dUj 

dx 


?) 


( 10 ) 


dJ\ _ 1 V 
dr 15 i/A 2 

In view of the form of the Rotta model, it is reasonable to seek terms at the next 
order proportional to UiUj. As in Sect. I, such terms do appear, but cancel exactly. 
This apparently ubiquitous cancellation was also obtained by Smith and Reynolds 18 
in an analysis of the s transport equation. Accordingly, the second order analysis 
in r] produces quadratic terms in the velocity gradients. Finite truncation of this 
series violates the requirement that return to isotropy be independent of the mean 
flow. Therefore, we must seek a reasonable approximate summation. The form of 
the lowest order term given in Eqs. (10) and (11) suggests 





Ti. 


Despite its triviality, this replacement does produce an approximate sum which 
agrees exactly with perturbation theory to lowest order. It therefore can be consid- 
ered a type of Pade approximation. 

At the infinite Reynolds number asymptotic limit 


n'a = ( 12 ) 

where, in the Yakhot-Orszag theory, Cr == Vje ~ 1.6. Equation (12) is therefore 
simply the standard Rotta model with Rotta constant ~ 1.6 in agreement with an 
earlier proposal of Yakhot 17 . 

A preliminary discussion of higher order summation may be appropriate. By 
analyzing the spectral dynamics of the return to isotropy, Weinstock 3 concluded 
that the shear and normal stresses relax at different rates. Although this behavior 
is obviously not accommodated by the Rotta model, it is consistent with the present 
theory: the perturbation series for IJ' is obtained from the series for r by multiplying 
the term of order n by the factor C n e/K for some constant C n . The C„ are all 
unequal; therefore, the Rotta model is not exact. Now comparison with the series 
for r shows 5 that relaxation of the shear stress is governed by the linear term S^, 
whereas relaxation of the normal stresses is governed by the quadratic term S^ 2 *. 
Since C 2 # C\ , these stresses relax at different rates. The difference is suppressed 
in the Rotta model, which arose in the present formalism by replacing all of the C n ■ 
by Cl 


2.3 Algebraic Reynolds Stress Models 

The approximation, due to Rodi 20 , of the Reynolds stress transport equation 
by an algebraic model under the conditions of semi-homogeneous flow (negligible 
diffusion of r and t/K approximately constant) takes the form 
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P-e 

K 


UiUj (0) = 


.( 0 ) 


dUj dUi\ tji , tj . 

dZ + p dx J +/?W V 


(13) 


where J7 and depend on t and V(7. Explicit solutions for r can be obtained, 
at least in principle, for any such approximation 25 . Briefly, one introduces a basis 
for polynomials in W, and VU T . The basis contains 11 terms of homogeneity 
order n < 5. Writing r as a sum of these terms with unknown coefficients and 
substituting in Eq. (13) leads to the explicit expression 

t/K = £ H i m) S i n) ( W ’ Vt/T ) (14) 

where H is a scalar function of VU and VU T such that 


H\ m) ~ |Vl7| m 

when IVl/l — ♦ oo. The assumptions made on the approximate summations require 
m + n = 0; thus, t/K is bounded when SK/e oo. For example, the familiar eddy 
viscosity formula is replaced in Eq. (29) by a term 


T ~ 


K 2 


(VC/,VC/ T ) (W + VZ7 r ) 


Pope observed 25 that the coefficients n) in Eq. (14) would certainly be in- 
tractably complex; although they could be expUcitly exhibited by symbolic com- 
putation, the result would only pertain to the particular implicit equation for the 
Reynolds stresses assumed initially in Eq. (13). Therefore, it is equally reason- 
able just to postulate simple forms for the functions This type of modeling 

could be particularly interesting when applied to the coefficients of the quadratica y 
nonlinear models of Refs. 5, 11, and 12. 


2.4 Discussion 

The present analysis of the Reynolds stress transport equation, based on the 
Yakhot-Orszag renormalization group and (tensorial) ^-expansion summation as 
suggested by Yakhot et al. 9 , has led to a model transport equation incorporating the 
well-known LRR and Rotta models. The analysis gives theoretical support both to 
these models and to the constants sometimes used with them. More significantly, it 
exhibits the LRR and Rotta models as lowest order approximations, and therefore 
also supports their replacement with higher order nonlinear models which would 
be deduced by more accurate approximate summations. The consistency of the 
analysis with higher order effects like the unequal relaxation rates of shear and 
normal stresses has been discussed. 


3. Future Plans 

The nonlinear eddy viscosity representation of the Reynolds stresses 
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2 K>,8Ui.au i \ 

r « = - a T + aS"J 



r Crl mv\c4^ + ^) w 

L \9xp dx p / Vaxp dxj ox p dxi / 



0Pp dUp 

dxi dxj 



in which if is the turbulence kinetic energy, e is the dissipation rate, C u , C r \ , C r2 , C t3 
are constants, and (0) denotes deviatoric part, was proposed by Yoshizawa 11 in or- 
der to model normal stress effects in shear flows by means of an explicit formula 
for the stresses. The significance of this formula is not limited to this property: 
Yoshizawa’s derivation using a special perturbation expansion, the two-scale direct 
interaction approximation, showed that the expansion could be continued to any 
order in the mean velocity gradient and thereby exhibited the Reynolds stress ten- 
sor as the result of an infinite number of increasingly complex interactions between 
the mean velocity field and turbulence. Related expansions are given in Refs. 5, 12. 

The infinite expansion which contains Eq. (15) can be written symbolically as 


r = KAqS° + —AiS 1 + ~ Y A * S i + 

6 £ i<N 3 


t<JV„ 


(16) 


Eq. (16) can be considered a decomposition of the Reynolds stress 

r = t° + r 1 -I 


(17) 


where 


T* = E = E T ” ( 184 > 

i<N n i<N n 

where S" denotes a symmetric tensor homogeneous of degree n in the mean velocity 
gradient VU and its transpose, N„ denotes the number of linearly independent 
terms of order n (so in Eq. (15), N 2 = 3), and the A tn are constants. 

In an analysis of the Reynolds stress transport equation by renormalization group 
techniques 35 , we found analogous expansions for the term which governs the return 
to isotropy, 


n' = ^[* 2 
k L £ 


£ £ z J 


i<N 2 


and for the rapid pressure-strain term 


(19) 


n = -k s 
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..£{c[ s . 


(0) 

vt^ + vt/s 1 ! + D\ 
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+^T £ [<** ( 5 < VpT + VUS?) m + D a (S?VU + Vl/ T S?)] + • • ■ (20) 

£ i<N 2 1 

Suppose that these sums are introduced into the Reynolds stress transport equation 

t = II' +U-P + D ( 21 ) 

where the dot denotes convective derivative, P = tVU t +VUt is production and D 
denotes the diffusion term. In order to obtain a Reynolds stress transport equation, 
it is necessary to express the sums (19) and (20) in terms of t and yU- t oug 
the coefficients A t B, C, D can be explicitly exhibited to any order in perturbation 
theory, they do not have an obvious law of formation. Therefore, the sums (19) 
and (20) can only be approximated by polynomials in r and VU if some hypotheses 
relating the coefficients is introduced. There is no unique hypothesis of this sort, 
but the simplest 35 seems to be 

Bi n /Bi = Ain/Au n> 2 (22) 

Cin/Ci = Din/Di = A in /A u n> 2 

which leads to the Rotta return to isotropy model and to an LRR model for the 
rapid term. The approximation expressed by Eq. (22) can be compared to the 
summation introduced in an analogous context by Yakhot et-al. 9 , and to Pad6 
approximation: it agrees with the perturbation theory of Eqs. (19), (20) to lowest 

order, but includes effects of all order in VU. 

By evaluating more terms of the perturbation series explicitly and introducing 
an equation like (8) for coefficients of higher order, a hierarchy of models could be 
generated. They would initially be nonlinear in VU , as advocated by Speziale , 
but one might introduce the perturbation series for r-r to obtain a model nonlinear 
in T 1,2 . However, the close analogy between Eqs. (19) and (20) and Yoshizawas 
expansion (16) suggests a different approach: namely, use Eqs. (17) and (18) to 
replace S t " in Eqs. (19) and (20) by r?/A in . Substitute these modified expressions 
and Eq. (17) into the transport equation Eq. (21), treat r” as having order | VU | n , 
and separate the terms of like order in | VU | in the standard perturbation theoretic 
fashion. The result is that the terms r t n in the decomposition (17), (18) themselves 
satisfy coupled transport equations. 

For simplicity, let us write an approximate system for the r n instead of for the 
r" and assume the most elementary scalar diffusion model. Then the the single 
transport equation for r would be replaced by a system 

f n = -C%r n + Ci (r n_1 VU T + VUr n_1 ) (0) 

+C 2 n (r n_1 VU + VU T r n-1 ) (0) + V 2 r n ,n > 1 (23) 

Since r = £ rT \ the system (6) should be constrained to satisfy Crow’s condition 
and to contain the exact production term following summation over n. Making the 
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coefficients C R , C”, C£, C£, independent of n reduces the system to a model of the 
LRR form for r = X) r n . 

This conclusion also follows from Leslie’s analysis 34 of the direct interaction (DIA) 
equations for shear flow. Leslie suggested a perturbative solution for the (tensor) 
correlation function and Green’s function 


U = U° + U l + 
G = G° + G 1 + 


(24) 


where U n and G n are of the order | VU | n , and observed that this expansion is 
simultaneously an expansion in powers of the mean strain, and a decomposition into 
symmetry types of increasing complexity. This is also a feature of the expansion 
(16). Substitution of Eq. (20) into the equations of the direct interaction approx- 
imation gives a coupled system for the U n and G n in standard fashion. Then in 
principle, by integrating each equation of this system over all wavenumbers and 
introducing the definitions 


-/ 


U n d k 


we could attempt to obtain coupled transport equations for the r n . Unfortunately, 
the derivation of equations for single-point quantities from DIA is not entirely 
straightforward, and more heuristic methods like two-scale DIA 8,14 and renormal- 
ization group are required. 

The simplest model of this type is a two component model, 


T 


-KI = t'+t 


2 


with 


f 1 = - CW - —K (VU + VU T ) + (Ci 1 - 1) (r 2 VU T + V*7r 2 ) (0) 
1 o 

+ C\ {t 2 VU + VU t t 2 ) {0} + C 1 d V^-Vt 1 
t 2 = - C 2 r t 2 + (C 2 - l) (t 1 VU t + VUt 1 ) (0) 

+ Cl { t 1 VU + VU 7 V 1 ) (0) + C 2 d V^-Vt 2 


The appearance of r 2 in the equation for r 1 is required by the production term 
constraint which is violated by direct truncation of the system (20) at the second 
order. The properties of this model for simple shear flow, in which dUi/dxj = 
S6n6ji follow by setting 


T 1 = 712 

1 

1 

T 2 = 

T U 

722 


0 


733 _ 


(26) 
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Then Cl describes the return to isotropy of the shear stress, and C\ the return 
to isotropy of the normal stresses. These relaxation rates can be unequal m this 
theory, an effect predicted by Weinstock’s analysis 3 of the spectral dynamics of the 
return to isotropy. Weinstock suggested further that the individual normal stresses 
relax at different rates: this effect is not accommodated by the present model, but 
could occur in a model in which r 2 is divided into the three tensor components 

Tl The inequality of the Rotta constants in the present model can be used to over- 
come a defect of models of the LRR type, that in semi- homogeneous flows m which 
a .. — t^/K is approximately constant, the ratio an is too small whereas an is 
too big 1 ' 0 . Following Speziale 19 , write Eq. (21) as a system of equations for the 
ratios ay and set a y = 0. There results, 


2anV = (CJi - 1) - j { C R ~ ^ [( C i ” l ) 022 + C * an " 15 1 } 


1/2 


an = 


u 22 = 


[-| (C 2 - 1) + §C 2 2 ]t/a 12 

+ 1) + t/ a l2 

[f (c? - 1) - 

(-C^ + 1) + t?ai2 


(27) 


When P/e = —V a n = 1> 


C\ — 1 C% 

a 22 ( 
'R U R 


fll2_ 15 C l R Cl an Cl au 


= [-| (Cf - 1) + |c|] /CS (») 

C22 =[| (C? - 1) - |cf]/ca 

Eq (28) shows that setting C\ < both increases a n and decreases a 12 and 
thus improves the agreement between theory and experiment. The trend required 
here is consistent with Weinstock’s findings 3 that the Rotta constant for the normal 
stresses should be smaller than the shear constant and should take values close to 

The inequality of the coefficients describing the rapid terms affects the behavior 
of the model under rapid distortion. The rapid distortion analysis of passive y 
strained turbulence predicts that a one-component limit state is reached in which 
r, = 2K6ii6n, V = SKje - oo and P/e is finite 21 . Whether or not this limit 
can in fact occur as an asymptotic state, a stress model should accommodate it 
because it can exist approximately in steady flows as a ’’spatial transient in strong y 
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inhomogeneous regions of very high convective or diffusive transport 17 . Moreover, 
models of the LRR form do not capture the transient evolution of rapidly distorted 
flows well. Let us assume that all model coefficients are functions of the ’’state” 
of turbulence, following Shih and Lumley 2 ; the precise parametrization of the state 
will be left to future investigations. From Eq. (13), it is evident that Pje = -rja^ 
will be of order 77 in the one-component state in which an = 4/3, an = -2/3 unless 


-I W - x > + H = is 


in this state. But Eq. (13) also shows that an = 4/3, a^ = -2/3 requires 


C\ ~ 0, C\ ~ 0 

in this state. These conditions are inconsistent in a model in which C{ = Cl and 
C\ = Cf> but are clearly consistent with the present proposal. 

The advantages of this model are consistency with a systematic perturbation 
theory, the possibility of unequal relaxation rates for normal and shear stresses 
in relaxing strained turbulence, improved agreement with experimental data for 
universal ratios in simple shear flows, and the possibility of accommodating the 
one-component limit of rapid distortion theory. 
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